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, Abstract 

The ground state configurations of the one-dimensional Falicov-Kimball model 
are studied exactly with numerical calculations revealing unexpected effects for small 
interaction strength. In neutral systems we observe molecular formation, phase 
separation and changes in the conducting properties; while in non-neutral systems 
^ . the phase diagram exhibits Farey tree order (Aubry sequence) and a devil's staircase 

structure. Conjectures are presented for the boundary of the segregated domain and 
the general structure of the ground states. 
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1 Introduction 



The study of the general properties of fermion systems, such as metal-insulator transitions, 
crystal formation and transitions to mixed- valence states, on the basis of microscopic mod- 
els, is a domain of vigorous research. Among the models considered nowadays growing 
attention is being paid to a simple model proposed almost a quarter of a century ago by 
Falicov and Kimball 0. This model was put forward to describe metal-insulator transi- 
tions in transition-metal and rare-earth materials, where an analysis of the experimental 
data suggested that these transitions are of purely electronic origin. 

The first investigations of the Falicov-Kimball model, carried out by means of approx- 
imate methods, like mean field or coherent phase approximations [0 , led to contradictory 
results and reduced the general interest in the model. However in 1986 Kennedy and Lieb 
H and independently Brandt and Schmidt ||] obtained rigorous results showing the exis- 
tence of a phase transition at sufficiently low temperatures. They considered the simplest 
version of the model with only one band and spinless fermions, later called the spinless 
Falicov-Kimball model. Even this simplified version has numerous, physically sound, in- 
terpretations ||, ||. In one of them, related to crystal formation, the system is viewed as 
consisting of quantum electrons hopping over a lattice and classical ions occupying some 
of the lattice sites. The interesting case occurs when the lattice consists of two equivalent 
sublattices and electrons are allowed to hop only from a site belonging to one sublattice 
to a nearest neighbour that belongs to the other sublattice. The only interaction is the 
on-site attraction (or repulsion) between electrons and ions and the hard core repulsion 
between ions. 

If the densities are equal to 1/2, or equivalently, if the chemical potentials correspond 
to the hole-particle symmetry point (half filled band case), then it has been proved that 
in dimensions two and higher this simple interaction leads to an order-disorder transition 
as the temperature is varied |§, |§ . 

In the low temperature phase the ions arrange themselves into a periodic structure, 
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where they occupy only one of the sublattices. An important point to achieve this remark- 
able result is the fact that, in distinction to other lattice fermion models, a study of the 
Falicov-Kimball model can be reduced to a study of tight-binding Schrodinger equations 
in a variety of potentials. At each site these potentials assume only two values, say or 
— U, that correspond to the absence or presence of an ion at this site. Let us underline 
that contrary to band theory, where the potential in the Schrodinger equation is fixed, the 
potential is here a variable and the problem is to find that potential whose contribution 
to the partition function is dominant (the so called annealed problem). 

One consequence of this fact is that the model can be studied by comparing the 
ground state energies corresponding to certain classes of potentials (ion configurations), 
i.e. by constructing restricted phase diagrams for a finite or infinite system. Initially the 
method of restricted phase diagrams was applied to a 2-dimensional system in || and to 
a 1-dimensional system in [H]. 

In recent years numerous papers devoted to the spinless Falicov-Kimball model ap- 
peared. They dealt mainly with the following topics:-various generalizations of the first 
rigorous results (quoted above) |7| , [U , U , JToJ] , -investigations of the segregation principle 
(first formulated by Freericks and Falicov [|]) [|ll|], fll2"l , -investigations of the effec- 
tive interaction between the ions [14fl, ||15|| , fl6| , [17] , -attempts to find out the general 
structure of the ground state phase diagram of the one- dimensional system, in particular 
new candidates for low-temperature configurations of the ions @, |18[]. We shall refer to 
the papers quoted here in due course. 

This paper belongs to the last of the mentioned groups. It is organized as follows. In 
section 2 we provide the reader with the relevant informations, definitions and notations. 
Then in section 3 we describe the results obtained in the canonical ensemble; after that 
in section 4 the results obtained in the grand canonical ensemble. Finally in section 5 we 
summarize our observations, conclusions and conjectures. 
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2 Definitions and notations 

The object of this work is to study the ground state properties of the one-dimensional 
spinless Falicov-Kimball model. This model describes a system of itinerant spinless quan- 
tum fermions interacting with classical ions on a one-dimensional lattice A. Whenever 
the lattice has a finite number |A| of lattice sites we impose periodic boundary conditions. 
The Hamiltonian is 



H A = Y^ (a* x ax+i + a* x+x a x - Uw{x)o* x a^j , (1) 

where the creation and annihilation operators of itinerant spinless electrons and 

w(x) = or 1 is the number of ions at the lattice site x. The total number of electrons 
N e = J2 a*a x and the total number of ions (in the configuration w) Ni(w) = w(x) 
are the conserved quantities. The function w(x),x G A, is called the ion configuration. 
The electrons can hop only between nearest neighbour sites and the corresponding kinetic 
energy matrix elements are set to unity. Therefore in there is only one independent 
parameter expressed in these units: the electron-ion interaction U. 

For a given configuration w(x) the Hamiltonian ([]]) is the second quantized form of 
the one-particle operator whose matrix elements are 



h X y = txy - Uw(x)S xy , (2) 

where 

1 otherwise. 

In the following, when we speak of the spectrum of the configuration w we mean the 
spectrum of this one-particle Hamiltonian. The bands (gaps) in this spectrum are referred 
to as the bands (gaps) of the configuration w. The ground state energy E u (w; N e ) of the 
system (P corresponding to N e electrons and the ion configuration w is equal to the sum 
of the N e lowest energy levels of the single-particle Hamiltonian @. 
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Applying to ([!]) the unitary hole-particle transformation with respect to the ions we 



obtain 



E iu) (w*; N e ) = E { - u) (to; N e ) - UN, 



(4) 



where w*(x) = 1 — w(x), while with respect to the electrons we get 



(w; N e ) = E { -^ (w; |A| - N e ) - UN t (w). 



(5) 



Therefore we can reduce the range of N e and A, considered and arbitrarily fix the sign of 
U (in the sequel we drop the superscript U by E). We choose U > 0, i.e. the particles 
attract each other at the same lattice site. 

We study also the infinite systems in the thermodynamic limit. The limit is taken in 
such a way that the particle densities per site A e /|A|, Aj/|A| and the ground state energy 
density E(w; N e )/\A\ tend to the electron density p e , the ion density and the ground 
state energy density e(w;p e ). For the periodic ion configuration w the Green function 
of the tight binding Schrodinger equation and consequently the density of states n(w), 
p e and e(w; p e ) can be determined exactly (T^]. However the corresponding formulae 
contain the band edges of the spectrum of w, which are zeroes of the polynomial whose 
order coincides with the period of w. At this point it is necessary to perform numerical 
calculations. 

The problem we want to investigate is to find the ground state configurations (g.s.c), 
i.e. the zero temperature phase diagram. In the canonical formalism one is given the 
electron and ion densities (p e , pi) and the task is to find the configurations w that min- 
imize the energy E(w,N e ), Ni(w) = Aj, or the energy density e(w,p e ), Pi(w) = p{. In 
the following sections we shall construct "restricted phase diagrams", which amounts to 
minimizing the set of functions e(w,p e ), where w runs over some specified class of ad- 
missible configurations, for example the class of all periodic configurations with period 
smaller than some fixed value. Since the minimum of a set of convex functions is not 
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necessarily convex, one has to take the convex envelope of the minimum. For those values 
p e where the convex envelope does not coincide with the minimum of the e(w,p e ), the 
g.s.c. is then a mixture of two or more configurations. The configuration w of ions in the 
mixture iui&u^ is defined as follows: The finite lattice A is partitioned into two parts Ai 
and A2; the restriction of w to Ai and A2 is W\ and W2 with |Ai|/|A| = a and 



The thermodynamic limit is then taken keeping a fixed. In particular the mixture of the 



clump together. The mixture wk,— will play an important role in our studies of neutral 
systems (section 3). The density of states of the mixture w&— is 



Pi = a piiwx) + (1 - a) Pi(w 2 ). 



(6) 



full and empty configurations is the so called segregated configuration^ where all the ions 



n a {n) = a n(w; //) + (! — a)n(— ; p) , 



(7) 



where 




and n(w; /j,) is given in []19| |. Consequently the electron and ion densities are: 



p e = a p e (w) + (1 - a) p e (-) 



Pi = a p/q , 



(9) 



where 




(10) 



—00 



and p e {— ) is given by eq. (|14D . Similarly the energy density of this mixture is 



e = a e(w) + (1 — a) e(— ) 



(11) 



where 




(12) 



—00 
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and e(— ) is given by eq.(|T5]). 

In the grand canonical formalism one is given the electron and ion chemical potentials 
(fjL e , Hi) and the task is to find the configurations w that minimize the free energy density 



/ (W, He, Hi) = e (W; He) - He Pe {w, He) ~ Hi Pi( W ) ■ ( 13 ) 

Again we shall construct "restricted phase diagrams" which amounts to minimizing 
f(w;He,Hi) over some class of configurations. In this case there is no need to take the 
convex envelope since the free energy density (|T3| ) is a concave function of the chemical 
potentials. 

In the following sections we shall consider special classes of configurations. It is useful 
to define them here and to introduce an appropriate notation. There are two translation- 
ally invariant configurations: the full configuration, denoted by +, where all the sites of 
the lattice are occupied (pi = 1) and the empty configuration, denoted by — , with no ions 
(pi = 0). The unique band of the empty configuration extends from —2 to 2 (while the 
band of the full configuration is translated by —U). At the Fermi level hf the electron 
density is 



pe (-; a*f) = < 



if hf < -2, 

vr" 1 arc cos (-fj, F /2) if \hf\ < 2, (14) 



1 if hf > 2. 

and the corresponding ground state energy density is 



-7i-y4-/4 if \hf\<2, 
< (-://,.) = { ' (15) 

otherwise. 

A periodic ion configuration that does not belong to one of the two classes considered, 
with period q and p ions per unit cell (p« = p/q) is denoted by p/q. 

The unit cell of the crenel configuration {p/q} consists of p consecutive sites occupied 
by the ions while the remaining q — p sites are empty. The atomic most homogeneous 



configuration with pi = p/q is denoted by \p/q}. In this configuration the ions in the unit 
cell are distributed in such a manner that the distances between two consecutive ions are 
either d or d + 1 with d < q/p < d + 1. Furthermore the distribution of the distances d 
and d + 1 has to be most homogeneous ||12|| , pOf . More precisely, with p relatively prime 
to q, the position of the ions in the unit cell is given by kj solutions of the equations 
pkj = j modq, j = 0, 1, . . .,p- 1. 

Similarly we define the n-molecule most homogeneous configuration [p/q] n , n > 2, 
where pi = p/q and p is a multiple of n. This configuration is defined in the same way 
as [p/q], but replacing "ion" by "n-molecule" that consists of n consecutive ions. More 
precisely, with p/n relatively prime to q, the position of the ions in the unit cell is given 
by kj solutions of the equations (p/n)kj = j modq, j = 0, 1, . . . ,p—l. In particular [n/q] n 
is identical with {n/q}. In the following we denote the n-molecules by H n . 

The following example illustrates the above definitions: for p = 6 and q = 10 the 
unit cell for the crenel configuration {6/10} is [•••••• o o oo], for the atomic most 

homogeneous configuration [6/10] it is [• • o • o • • o«o], i.e. it is the extension of the unit 
cell corresponding to [3/5], for the 2-molecule most homogeneous configuration [6/ 10] 2 
it is [• • o o • • o • «o], for the 3-molecule most homogeneous configuration [6/10]3 it is 
[•••oo»»»oo]. Finally as an example of a periodic configuration 6/10 that is not a member 
of the above mentioned classes we can take the one with the unit cell [• • o • • o • o »o]. 
The list of all the periodic ion configurations that enter into our phase diagrams is given 
in Table ^ The mixture of the configuration w = p/q and the empty configuration is 
denoted p/qh,— or p / q&ivacuum (if the Fermi level pp is smaller than —2). 



S 



Pe 


Pi 


unit cell 


1/11 


[2/H]2 






[3/H] 3 




1/10 


[2/10] 2 






[3/10] 3 






[4/10] 4 






[5/10] 5 




1/9 


[2/9] 2 






[3/9] 3 






[4/9] 4 




1/8 


[2/8] 2 






[3/8] 3 






[4/8] 4 




1/7 


[2/7] 2 






[3/7] 3 






[4/7] 4 




1/6 


[2/6] 2 






[3/6] 3 




2/11 


[4/ll]a 


••ooo»»oooo 


1/5 


[1/5] 





Pe 


Pi 


unit cell 


1/5 


[2/5] 2 






[3/5] 3 




4/19 


[4/19] 




3/14 


[3/14] 




2/9 


[2/9] 


•oooo»ooo 




3/9 


••ooo»ooo 




[4/9] 2 


MOOMOOO 




5/9 


IIOOKtOO 




[6/9] 3 


• • • o • • • c o 


1/4 


[1/4] 






[2/4] 2 




3/11 


[3/H] 


•ooo«ooo«oo 




4/11 


••oo»ooo»oo 




5/11 


••oo»»oo»oo 




[6/H]2 


•IOOIIOOIIO 




7/11 


• • • O • • O • • o 


2/7 


[2/7] 


• O • 




3/7 


• • O • 




[4/7] 2 


• • O • • 



Table 1: (a) The periodic ion configurations that appear in the restricted phase diagrams. 
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Pe 


Pi 


unit csll 




Pe 


Pi 


11 Tilt Ppll 


2/7 


[5/7] 


• • • o • • o 




5/12 


8/12 


• • • • o • o • • o »o 


3/10 


[3/10] 


• OOOIOOIOO 






9/12 






4/10 


• IOOIOOIOO 




3/7 


2/7 


• O • O O o 




5/10 


• IOOMOIOO 






[3/7] 


• O • O • o 




[6/10] 2 


• • • • • • 






4/7 


• O • • • o 


1/3 


[1/3] 








5/7 


• O • • • • 




[2/3] 2 






4/9 


3/9 


•o»o»oooo 


3/8 


2/8 


• O O • o oo 






[4/9] 


•o«o«o«oo 




[3/8] 


• O O • o »0 






5/9 


•oaoaoaao 




4/8 


• • O • O »0 






6/9 


• O • O • • • • 




5/8 


• • • • o »0 




5/11 


3/11 


•o»o»oooooo 




[6/8] 2 








4/11 


•o»o»o»oooo 


2/5 


1/5 


• O O O 






[5/11] 


•oaoaoaoaoo 




[2/5] 


• O O • 






[6/11] 


•oaoaoaoao* 




3/5 


• • o • o 






7/11 


• • • o • o • o • o • 




[4/5] 2 


• • • • 






8/11 


• • • • • o • o • o • 


5/12 


[5/12] 


•oo»oo«o«o«o 




1/2 


[1/2] 






6/12 


• OOIOIIOIOIO 












[7/12] 


••o»o»»o»o»o 











Table 1: (b) The periodic ion configurations that appear in the restricted phase diagrams; 
continuation. 
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3 Neutral systems. Phase diagrams in the canonical 
ensemble 



In this section we study the ground state properties of the neutral systems in the canonical 
ensemble. If the system is finite we thus take N e = N{ = N and we are interested in the 
configurations w that minimize the energy E(w; N). 

If the system is infinite, then p e = Pi = P and the problem is to find the configura- 
tions w that minimize the energy density e(w;p/q). We shall first consider only periodic 
configurations with p = p/q, p and q being relatively prime integers. 

For large U and rational densities it was first conjectured H and then rigorously proved 



that the g.s.c. is the most homogeneous configuration of atoms [|I2 |. Since this result was 
established using the perturbation series of h{w) in powers of U~ x , the question naturally 
arises whether the restriction to large U is technical or fundamental. Until now it has 



been usually conjectured that such a theorem should be valid for all values of U |T2|, |2T 



Indeed, the results of Kennedy and Lieb || and Brandt and Schmidt M for the half-filled 
band case (p = 1/2) and the numerical studies of Freericks and Falicov || for pi = 1/2, 
1/3 seemed to support this conjecture. 

However the analysis of a finite number of ions and electrons on an infinite lattice 
(with N e = Ni), while providing a qualitative explanation of this result, suggests also 
that it cannot remain valid for small U. Indeed, for large U and to leading order in U^ 1 
the energy of the neutral system is given by an effective two-body potential of the form 
2(d+ 1) exp[— X(2d + 1)], where d is the distance between the ions fllEf . In other words, 
for large U an equal number of ions and electrons will form neutral atoms which repel 
each other with an effective two-body potential that is convex and decreasing. Using the 
result of Hubbard ptl we thus conclude that for large U the g.s.c. of the neutral system 
with a finite density should be the most homogeneous configuration. 

On the other hand, one can easily check that for U < 2/ the energy of the neutral 
2-molecule is smaller than the energy of two neutral atoms that are infinitely separated. 
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Indeed the ground state energy E atom of the neutral atom is 



E 



'atom 



(16) 



while the energy E 



'2— molecule 



of the neutral 2-molecule is 



U 3 



if U>2, 



-2 



U 2 - 1 
(U + 2f 
U + l 



E 



2— molecule 



= < 



(17) 



if U <2. 



Moreover for U < 2 the 2-molecule has only one bound state. 

More data are contained in Fig. 1, where the energy of two ions and two electrons on 
the infinite lattice is shown as a function of the distance between the ions for different 
values of U. For large values of U the ground state consists of two neutral atoms at infinite 
distance. On the other hand, for small values of U this happens when one 2-molecule 
with one bound electron and one diffusion electron is formed. Similarly, Fig. 2 shows that 
the ion configuration that minimizes the energy of three ions and three electrons depends 
on the value of U. There are three kinds of ground states: 

- three neutral atoms at infinite distance (U large), - one 2-molecule with one bound 
state and one neutral atom at infinite distance, - one 3-molecule with one bound state 
(U small). 

Finally we can see in Fig. 3 that the energy between two 2-molecules (and 4 electrons) 
presents the same properties as the energy between two ions: for sufficiently large U it 
is convex and decreasing and thus the ground state consists of two 2-molecules infinitely 
separated, while for U sufficiently small the ground state consists of one 4-molecule with 
one bound state. Although we have not checked it, we expect that the energy between 
two n-molecules (and 2n electrons) has similar properties. To conclude the above consid- 
erations of neutral systems with a finite number of ions, let us find the n-molecule with 
the lowest energy per ion as a function of U. The result is shown in Table |2|. We note 
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Figure 1: Energy per ion for two ions and two electrons versus the distance between the 
ions. 



Figure 2: Energy per ion for three ions and three electrons versus the two distances 
between the ions. 
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Figure 3: Energy per ion for two 2-molecules and four electrons versus the distance 
between the molecules. 



[/-interval 


Type of n-molecule 


(1.154, oo ) 


#i 


(0.511, 1.154) 


H 2 


(0.279, 0.511) 


H 3 


(0.175, 0.279) 


H 4 


(0.119, 0.175) 


H 5 


(0.086, 0.119) 


H 6 


(0.065, 0.086) 


H 7 


(0.051, 0.065) 


H 8 


(0.041, 0.051) 


H 9 



Table 2: n-molecules with the smallest energy per ion versus U. 
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type of molecule 


unit cell 


#1 


(.)(13o)(.)(13o)(.)(13o)(.)(14o)(.)(13o)(.)(13o) 
(.)(13o)(.)(13o)(.)(14o) 


H 2 


(2«)(26o)(2«)(27o)(2«)(26o)(2«)(27o)(2«)(26o) 
(2«)(27o)(2«)(26o)(2«)(27o)(2«)(26o) 




(3«)(39o)(3«)(40o)(3«)(40o) 


#4 


(4.)(52o)(4«)(53o)(4«)(53o)(4«)(53o)(4»)(53o) 
(4«)(53o)(4«)(53o)(4.)(53o)(4«)(53o) 



Table 3: Periodic configurations for p = 9/128 taken into account in Fig. 4. In our 
notation, for instance (3») stands for three consecutive sites occupied by ions while (3o) 
stands for three consecutive empty sites. 

that for each n-molecule there is an interval of U where its energy is minimal. Moreover 
in this interval the n-molecule has only one bound state. We then conjecture that for 
given U, the neutral system with zero density of particles will form n-molecules, with 
n = n(U), which are infinitely separated. 

The above results suggest that, for sufficiently small U, the neutral systems with 
nonzero density might realize their ground state by forming n-molecules, n > 1, that are 
most homogeneously distributed. To test this idea we have considered a finite system of 
512 sites and studied the restricted phase diagram for configurations of the form [p/q] n . 
For example, for pi = 9/128 (iVj = N e = 36) we have considered the ion configurations 
whose unit cells are displayed in Table [| The results are given in Fig. 4. To determine 
the finite size effect, the transition line between the atomic configurations Hi and the 2- 
molecule configurations H 2 has been computed for a system of 200 sites and no significant 
difference was observed. We also compared the energy of the finite system with the energy 
of the corresponding periodic configuration of the infinite system (e.g. see Table |5] and 
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Figure 4: Restricted phase diagram with respect to n-molecule most homogeneous con- 
figurations, n < 4, over the lattice of 512 sites. The points are special values discussed in 
the text. 

Table |||). From this analysis we can thus conclude that the usual conjecture concerning 
the ground state of a neutral system (mentioned at the beginning of this section) cannot 
be correct: for sufficiently small U and sufficiently small densities the most homogeneous 
configuration of atoms cannot be the ground state configuration. 

Of course from this analysis we cannot conclude that the molecular configurations 
[p/q]n are the true g.s.c. Therefore for some special values of U and p (Fig. 4), and for 
an infinite system, we have computed the energies of all the periodic configurations with 
period smaller or equal to some Q. The configurations with the smallest energy and 
their energy are given in Table 

Clearly, for sufficiently small U and p the most homogeneous configurations of atoms 
or molecules are not the g.s.c. These results show that theorem 1 of Freericks and Falicov 
|| is not always correct: for p = p/q there exist periodic configurations whose energy 
is smaller than the energy of atomic configurations and whose period is larger than q 
(similar results hold true for smaller values of U). However, we shall see in the sequel 
that, if the true g.s.c. is periodic then it is given by the construction of Freericks and 
Falicov || and this remains true in the non-neutral case (p e ^ Pi). 
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u 


p 


period < Q : 


unit cell 


ener£ 


;y 


energy (fl/pl) 


o in 


1/5 


20 


• o o o o 


- 378 


561 




n m 


1 /6 

l/U 


24 


• •nnnnii I 1finl 

W W 1 — > \ ' * ' W W W 1 _[_ \ 1 1 — I 1 


- 321 


41 3 


- 321 386 


10 


1 /7 


21 


• • •( 1 8o ) 


- 278 


539 


- 278 51 1 


0.20 


1/6 


24 


iiooooii (16o) 


- 0.325 


073 


- 324 972 


20 


1/7 


21 


• •ooooo© ^13o^ 


- 281 


342 


- 281 238 

U. ZU ± iuU 


50 


1/5 


20 

Zj W 


• o o o o 


- 401 


926 

Zj \j 




0.50 


1/6 


24 


• • o o o • •fl7o) 


- 0.338 


876 


-0.338 655 


0.50 


1/7 


21 


• •ooooo* (13o) 


- 0.292 


371 


-0.291 988 


0.54 


1/6 


18 


• • o o o o #(110) 


- 0.340 


847 


-0.340 815 


0.60 


1/6 


18 


• 00 


- 0.344 


206 




1.00 


1/6 


24 


• 00 


- 0.371 


472 




1.00 


1/7 


21 


• 


- 0.318 


390 





Table 4: The canonical restricted phase diagram for an infinite chain and special values of 
(U, p). The ground state configurations were selected among all the periodic configurations 
whose period does not exceed Q. 



u 


unit cell 


energy 


0.1 
0.2 


• • (5o) ••(5o) ••(26o) 

• • (4o) ••(4o) ••(28o) 


- 0.278 265 

- 0.281 258 



Table 5: Ground state configurations of the neutral system with 42 sites and density 
p = l/7. 



Furthermore Table |4] shows that for small densities and small U the particles tend to 
arrange themselves into molecules (as it was the case at zero density) and lead us to think 
that the g.s.c. could be mixtures of periodic configurations of molecules and the empty 
configuration. To check this guess we have considered a system of 6 ions and 6 electrons 
on a lattice of 42 sites, i.e. p = 1/7. The energies of all ion configurations have been 
computed and the g.s.c. selected (see Table We see that the g.s.c. of this finite system 
(p = 1/7) are the mixtures of {2/7}&— or {2/6}&— . However Fig. 4 suggest that the 
mixture of {3/6}&— could come into competition if the lattice were sufficiently large. We 
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u 


configuration 


energy (L = 420) 


energy (L = oo) 


0.1 


{2/7} & - 


- 0.278 580 


- 0.278 679 




{3/6} & _ 


- 0.278 538 


- 0.278 535 




{2/6} & - 


- 0.278 500 


- 0.278 500 




[1/7] 




- 0.278 511 


0.2 


{3/6} & - 


- 0.281 470 


- 0.281 490 




{2/6} & - 


- 0.281 452 


- 0.281 455 




{2/7} & - 


- 0.281 441 


- 0.281 439 




[1/7] 




- 0.281 238 



Table 6: Energy density of the mixture {p/g}&— for the neutral system with density 
p = 1/7, for a lattice of 420 sites and for an infinite lattice, compared with the energy 
density of the atomic most homogeneous configuration [1/7]. 

have thus computed the energies of these three mixtures for a lattice of 420 sites and for 
the infinite lattice (see Table |6|). These new and surprising results show that for sufficiently 
small values of U and p the g.s.c. is indeed a mixture of some periodic configuration with 
the empty configuration. We have thus looked at the restricted phase diagram for all 
configurations that are mixtures of the empty configuration and a periodic configuration 
with period q < 7 for p = 1/5,1/6,1/7 and period q < 10 for p = 1/10. Note that, 
in particular, the most homogeneous configurations are among those considered. Table 
[7| shows that the configurations minimizing the energy are mixtures of {p/q}&— (i.e. 
crenel configurations followed by the empty configurations) or the most homogeneous 
configurations of atoms. This observation suggests to restrict our attention to those 
configurations that are either most homogeneous or mixtures of the form {p/q}&— ■ We 
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p\u 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1/5 


{2/5} 


{2/5} 


















#i 


1/6 


(2/61 
1^7 u x 


(2/61 


12/5} 


12/5} 


12/5} 


12/5} 


12/5} 


Ha 

±J - 1 


Ha 


Ha 


JJ 1 


1/7 


{2/7} 


{3/6} 


{2/6} 


{2/6} 


{2/6} 


{2/6} 


{2/5} 


{2/5} 


{2/5} 


#1 


#1 


1/10 


{3/9} 


{3/8} 


{3/8} 


{2/8} 


{2/7} 


{2/7} 


{2/7} 


{2/6} 


{2/6} 


{2/6} 


#1 



Table 7: Restricted phase diagram for mixtures of the form p/qSz— (infinite lattice). 



have thus constructed the restricted phase diagram for all configurations that are either 
most homogeneous or mixtures of {p/q}Sz— with q < 10. This has been done in the 
range U G [0.01, 1.2], p G [0.1,0.3], taking steps AU = 0.05, Ap = 0.005. We have also 
considered the values of p (in the above range) of the form a/b with b < 29. In the domains 
that have not been expected to appear (e.g. {3/6}, between {2/6} and {2/7}) finer steps 
were used: AU = 0.01, Ap = 0.001. Finally the transition lines separating different 
domains of 2-molecules have been computed explicitly. The results are summarized in 
Fig. 5. 

To achieve a better understanding of the transition we have fixed U = 0.6 and we 
have found that for p G [0.14,0.15] the mixture of {2/5}&{2/6}&— has lower energy 
than {2/5}&- or {2/6}&-. Similarly for p G [0.18,0.19] the mixture {l/5}&{2/5}&- 
has smaller energy than {l/5}&— or {2/5}&— . Therefore the transition lines should 
be replaced by transition boundaries with mixtures of one or two periodic molecular 
configurations and the empty configuration. Our investigations of the grand canonical 
phase diagram, described in the sequel, suggest that between two domains \p/q] n Sz— and 
\p'/q'] n Sz— there is the domain [{jp +p')/(q + q')] n Sz— and so on. On the other hand 
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Figure 5: Restricted phase diagram of the infinite neutral system in the canonical en- 
semble. All the mixtures of the form {p/q}&— , with q < 10, and the atomic most 
homogeneous configurations whose period does not exceed 10 were taken into account. 
The dots correspond to the entries of Table 7. The configurations found in the grand 
canonical analysis are located at the vertical segments. The fat vertical segments repre- 
sent the mixtures {l/5}&{2/5}&- and {2/5}&{2/6}&-. The domains labeled 1 and 2 
are [4/9] 4 and [4/8] 4 , respectively. 

between two domains \p/q] n &— and \p' / q'] n+1 k,— there is a mixture [p / q] n &,[p' / q'] n+ i&L— . 

Let us summarize our observations. For neutral systems (p e = pi = p) and for U < 
2/\/3 there is a critical density p c (U) such that for p > p c (U) the g.s.c. is atomic most 
homogeneous (as for U > 2/y/3) and thus the g.s.c. varies from point to point. On 
the other hand for p < p c {U) larger and larger molecules appear, the g.s.c. remains a 
mixture of the periodic n-molecule configuration [p/q] n and the empty configuration, over 
a certain range of densities. As p decreases the period q increases discontinuously until a 
value where a mixture [p/q}( n +i)&— takes over. 

Furthermore, if the g.s.c. is atomic most homogeneous, then the system is an insulator, 
while if the g.s.c. is a mixture, the system is a conductor. We thus have a transition from 
an insulator to a conductor as U decreases. 

We now give an interpretation of our observations. We have already noticed that a 
finite number of ions and electrons on an infinite lattice will tend to make n-molecules 
(with one bound state only) which repel each other. For finite densities Pi = p e = P 
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Figure 6: Band structure of the configurations {2/5} and {2/6} for U between 0.1 and 
1.0. The vertical lines show the location of the Fermi level for the densities p = 1/10, 1/7 
and 1/6. 

we would like to understand why the g.s.c. is a mixture of {n/g}&— and what is the 
physical mechanism responsible for the period q. Looking at the band structures of the 
periodic configurations {n/q} we notice that the Fermi level p F of the mixture {n/g}&— , 
corresponding to p, is in the gap between the first two bands (Fig. 6) and thus we have 

Pi = p = a-, p e = p = a- + (1 - a)p e (-; p F ) • (18) 
q q 

Thus 

= a± = £, p e (-;p F )=p^-, (19) 
q n n — pq 

i.e. there is one electron per molecule in a bound state associated with the periodic 
ion configuration while the remaining p{n — l)/n electrons are in the extended states 
associated with the empty configuration of the ions. The corresponding energy density of 
the system is 
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e =£[g C W( {n / ff })_2(n-l)M , (20) 

" V y J 

where eP^({n/q}) is the energy density of the first band of the configuration {n/q} and 
y = irp(n — l)/(n — pq). For a given p and n the first term decreases as q increases, 
which can be seen as a consequence of the effective repulsion between n-molecules due 
to the bound electrons. On the other hand the second term increases as q increases, 
which reflects the fact that the electrons in the empty space tend to increase the volume 
occupied. Therefore the period q that minimizes the energy density results from the 
competition between the effective repulsive force between the molecules and the pressure 
exerted by the electrons in the empty configuration of the ions. 

4 Phase diagrams in the grand canonical ensemble 

To obtain a better understanding of the phases at zero temperature, and to avoid the 
difficulties associated with mixtures that one has to consider in the canonical approach, 
we shall now investigate the infinite system in the grand canonical formalism. The problem 
is the following: given the chemical potentials (p e ,pi), what is the configuration w that 
minimizes the free energy density f(w;p e ,pi). 

In the following we consider only periodic configurations and we construct the phase 
diagrams restricted to periodic configurations with period smaller or equal to 10, for inter- 
action U taking values 0.6, 0.4, 0.2. Using the symmetry properties and the known results 



15| we can restrict our investigation to the region p e 6 [—2 — U, —U/2], Pi G [—£7/2,0]. 
We have thus constructed the restricted phase diagrams in the following manner: For 
fixed p e , and given w, the free energy is a linear function of p^. We can therefore calculate 
exactly the interval of pi for which the configuration w yields the minimum of the free 
energy (among all the configurations whose period does not exceed 10). This was done 
taking steps Ap e = 0.05. Near the end point (p e = pi) and near the transition lines finer 
steps were taken. The resulting phase diagrams are presented in Fig. 7. To study these 
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diagrams we have also computed the band edges of all configurations that appear in Fig. 
7 and we have marked the intervals of \x e in which they are g.s.c. (see Fig. 8). 

Our analysis of the restricted phase diagrams, together with the information on the 
band structure, shows that the (p e , pi) plane is divided into three main domains: D + , 
where the full configuration is the g.s.c, _D_, where the empty configuration is the g.s.c. 
and D, where the g.s.c. is different from the full or empty configuration. In the range 



/i e < —U — 2 we already know |L5| that for /ij > the g.s.c. is the full configuration, for 
fii < the g.s.c. is the empty configuration and for /z; = all configurations have the 
same energy. For —U — 2 < ji e < /!*, where p* e < —2 (see Fig. 7), the domains D + and 
D- are separated by the curve 



f(+;H e ,Hi) = f[-]fJL e ,m) = 0, (21) 



i.e. 



Tc/Ji = — W A — nl — n e arc cos (— jl e /2) , 

(22) 

A point (fjL e , Hi) on this curve corresponds to a mixture of the full configuration with 
electron density 



p e (+,/i e ) = 7r 1 arc cos (-/i e /2) (23) 

and the empty configuration with p e (— ; /i e ) = (vacuum). This is the so called segregated 
configuration, in which all ions clump together ||. Along the curve (p^) we have 

Pe/ Pi < Pe ( + ; H* e ) ■ (24) 

In other words for densities p e and pi satisfying inequality (p^ ) the g.s.c. is the segregated 
configuration. The problem of determining the value of p* e will be addressed later on. 
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Figure 7: Restricted phase diagram of the infinite chain in the grand canonical ensemble. 
All the periodic configurations whose period does not exceed 10 were taken into account. 



24 



Figure 8: Bands of the g.s.c. that appear in the grand canonical phase diagram (Fig. 7) 
for U = 0.6. The continuous horizontal lines mark the extents of the bands, the dotted 
segments the /i e -extents of the corresponding g.s.c. 
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Fig. 8 shows that for each configuration w appearing in the phase diagrams (Fig. 7) 
the Fermi level is located in a gap of the corresponding spectrum. This property has 
also been observed for U = 0.4 while for U = 0.2 it does not seem to be true since the 
boundary between two domains with p e — p/q and pi — p'/q, ij>' + l)/q is not a straight 
line. However we have checked that this property is indeed verified when we consider 
configurations with larger period. 

Fig. 7 shows that the domain D consists of connected subdomains D( PejP .) charac- 
terized by some periodic configuration of the ions and by a fixed density of electrons p e 
(since the Fermi level is in the gap). We note that for given densities (p e ,pi), if there is 
a periodic g.s.c. corresponding to this density, then it is associated with one connected 
domain .D( PejP .). This is a consequence of the fact the the free energy surface is concave 
and at zero temperature it is piecewise afline. It turns out also that the configuration of 
the ions can be determined by the procedure given by Freericks and Falicov ||: 

Let p e = p/ q, with p relatively prime to q and pi = Pi/q, where Pi may not be relatively 
prime to q. Define the integer kj by the relation pkj = j mod q,j = 0, 1, 2, . . .. Then the 
position of the ions in the unit cell is given by kj with j = 0,1,2, ... ,pi — 1. 

The following structure of D is remarkable: the union D Pe of the domains D^ PetP .^ 
corresponding to the same electron density p e = p/q constitutes a connected region. In 
this region the ion density takes on the values p'/q, (p' + l)/q, . . . ,p"/q, for some p' and 
p" > p' (see Fig. 9). Moreover the unit cell of the configuration associated with the 
ion density pi = {p + l)/q can be obtained by adding one ion to the unit cell of the 
configuration associated with pi = p/q. 

Another property of D is related with the observation that the domains D n , where 
the ion configurations consist of n-molecules with fixed n are connected regions whose 
size decreases as n increases. There is a relatively large domain D\ corresponding to 
neutral atoms (p e = pi = p) with p ranging from some p m in(U) to 1 — p m in(U), a smaller 
domain D 2 corresponding to 2-molecules with pi = 2p e , and consecutively domains D n 
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Figure 9: Domains D Pe with p e = 3/8, 2/5, 3/7 and 4/9. 

corresponding to n-molecules, where n = 2,3, . . . ,N(U), with = np e (i.e. in all the 
domains D n there is one electron per molecule). In each domain D m ,m = 1, . . . ,N(U) 
the ion density satisfies the inequality p^i n (U) < pi < p^ x {U). Furthermore the minimal 
value Pmin{U) is zero for all m between some n(U) and N(U) (in the cases considered 
we had n(U) = N(U) — 1) and for those m the boundary between D m and the vacuum 
(p e = pi = 0, i.e. p e < —2) is linear. Between two domains D n and D n+ i the periodic 
configurations are such that their unit cells contain n and (n+l)-molecules only. Between 
D n and the vacuum the unit cells consist of molecules and a sequence of empty sites while 
between D n and the full configuration the unit cell contains atoms and a sequence where 
all sites are occupied. 

We note also that all the domains D Pe (with a given electron density), that appear in 
the restricted phase diagrams, follow the Aubry sequence (Farey's tree, Fig. 10) |22| with 
respect to p e and exhibit a devil's staircase structure with respect to their size. 

We therefore arrive at the conjecture that to all electron densities between zero and 
one there correspond, according to the Aubry sequence, the domains D Pe which present 
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Figure 10: Farey tree (or Aubry sequence) for electrons. 



a fractal structure known as the devil's staircase ]23 . 

Similar observation can be made with respect to the structure of the domains D n of 
n-molecules. The ion densities in these domains appear to follow the "n-molecule Aubry 
sequence", that is similar to the Aubry sequence, but starts with 0/n and n/n and one has 
to consider rationals p/q with p = mod q (see Fig. 10). This property can be understood 
as resulting from an effective two-body interaction between n-molecules , that is convex 
and decreasing (as was the case for the effective two-body interaction between the atoms 
for large U). However, because of the effective fc-body interaction, we do not have the 



exact devil's staircase of Burkov |p2 |. In the neutral case this effect was already observed 



in 



18| , where the periodic g.s.c. up to period 6 and for large U were studied. To check 
the above conjecture we have taken into account some periodic configurations with period 
larger than 10. The results are shown in Fig. 11. 

An inspection of the phase diagrams (and Fig. 12) yields a more detailed picture of 
the neutral case discussed within the canonical formalism. For U < 2/\/3 the domain 
corresponding to non-trivial neutral configurations (p e = consists of the curve AB, 
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Figure 11: Test of the conjecture regarding the Aubry sequence and the devil's staircase, 
(a) U = 0.6, (b) U = 0.6, (c) U = 0.2, (d) /z e -extents of the domains contained in D x . 
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part of the boundary between the empty configuration and the periodic configurations, 
and the domain Di ending at B. Along the curve AB the density increases monotonically 
from zero to PminiU). At each point on these curves the g.s.c. is a mixture \p/q] n &— or 
[p/g] n &[p'/g'] n+ i&— . Furthermore at each point the density is different (because of the 
weight of the vacuum, not because of the periodic g.s.c). In the domain D\ the density 
is piecewise constant, increasing from p m i n (U) to 1 — p m in(U), with a devil's staircase 
structure. The results for p e = pi and U = 0.6, 0.4, 0.2 have been reported in Fig. 5. 

We conclude this section with an argument leading to an equation for /!*, the p e - 
coordinate of the point, where the coexistence line of the empty and full configurations 
joins the domain D. The reduced phase diagrams show that for p e e [pi, —2] the bound- 
ary of D with the vacuum (p e — pi — 0) is piecewise linear and corresponds to the 
configurations of the form [n/oo] n . On the other hand our considerations in the canonical 
formalism have shown that for a finite number of ions the g.s.c. of a neutral system con- 
sists of n-molecules (where n is U-dependent) with one bound electron and n — 1 diffusion 
electrons. For p e < —2 there is no diffusion electrons and we then have only n-molecules 
with one electron. The free energy per ion is thus 

pW _ 

f{H n -p e) p t ) = ^ ££-AH. (25) 

n 

where Eq 1 ^ is the energy of the n-molecule ground state. On the other hand the free energy 
of the vacuum is zero. Therefore the boundary of the vacuum domain (for p e < —2) is 
given by the concave envelope of the functions of p e 

f(+;p e ,Pi) = 0, /(#„; p e , pi) = 0, i.e. p™ = (E^ — p e ) /n. (26) 

For p e = —2 the value n(U) that corresponds to the minimum of the sequence (Eq 1 ^ + 2)/n 
is equal to the number of ions in the molecule at zero density (Table 0). Since the function 
given by f(H n ;p e ,pi) = has the slope — 1/n, it is sufficient to construct the concave 
envelope of the functions given by eq. (|26|) with n > n(U). These functions are plotted 
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Figure 12: The region of the (/x e , /li) plane where the system is neutral (p e = pi). In 
the hatched area the g.s.c. are periodic while at the curve AB the g.s.c. are mixtures of 
periodic configurations with the empty configuration. 
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Figure 13: Graphs of the functions defined by eq. (26). 

in Fig. 13. In particular it appears that for U > 2.65 the value of \x\ is determined by 
f(Hi;/i e ,iii) = 0, i.e. 

fjLi = -V£+lP-n e , (27) 
which together with eq.(|22|) gives the following equation for /i*: 



(28) 



V4 + U 2 = U+^ (Jl - Jl - %arc cos (fi e /2)\ 



jl e = /j,* e + U, (*7>2.65). 
and the range of densities corresponding to the segregated configuration: 



Pe < Pi-arc cos (-p e /2) . (29) 

7T 

This curve is plotted in Fig. 14 together with the upper and lower bounds derived by 



Brandt [13]. To the precision of Fig. 14 the curve obtained by Freericks and Falicov 



almost coincides with our curve and only a few points have been reported. 
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Figure 14: Canonical phase diagram. For U > 2/\/3 and p e = pi the g.s.c. is the atomic 
most homogeneous configuration. The boundary of the segregated configuration, shown 
as the continuous line, was computed from eq. (29). The dashed curves represent the 
bounds obtained in [13] while the crosses represent numerical results of [6]. 

Finally we add a remark concerning the canonical phase diagrams obtained by Fre- 
ericks and Falicov 0. Let us take U = 0.6. For the densities (p e = 1/3, pi = 1/2) the 
unit cell of the g.s.c. in |J consists of three occupied sites followed by three empty sites; 
the energy density of the corresponding g.s.c. is equal to —0.667401. We found that the 
even mixture [l/3]&[2/3] has the energy density —0.668167. Similarly, for the densities 
(p e = 1/2, pi = 1/3)) the unit cell of the g.s.c. in |j] has the length 6 and consists of 
two occupied sites separated by one empty sites, the other sites being empty. The energy 
density of the corresponding g.s.c. is equal to —0.753835. We found that the mixture 
[l/2]&— has the energy density —0.754651. 

In conclusion, to be physically meaningful, an incoherent phase diagram (in the plane 
(p e , U) and for fixed p, ) must take into account mixtures of configurations with different 
ion densities. The same remark, as we have already seen, applies to the neutral case. 
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5 Conclusions 



Our analysis of the zero temperature restricted phase diagrams of the spinless Falicov- 
Kimball model has revealed several new and surprising effects. 

In the case of neutral systems we have seen that for U > 2/y/3 the ground state 
configuration is, as expected, atomic most homogeneous and the system is an insulator. 
However for U < 2/y/3 this remains true only for densities p which are not too small 
or not too large. Indeed there is a critical density p c (U) such that for p < p c (U) (or 
p > 1 — p c (U)) the ground state presents entirely different properties. We observed three 
new phenomena: 

-Phase separation; the ground state configuration is a mixture of a periodic config- 
uration and the empty configuration, i.e. all ions are distributed periodically in one 
semi-infinite half-space while the other half-space is free of ions. 

-Formation of n-molecules; the ions form clusters where n consecutive lattice sites 
are occupied and these clusters are most homogeneously distributed over the occupied 
half-space. 

-Changes in conducting properties] the Fermi level is located in the band of the empty 
configuration and therefore an extra electron can be added at no cost in the energy. 

As the density p decreases the following infinite sequence of phase diagram transfor- 
mations is observed. In the first stage the weight of the half-space with no ions increases 
while the unit cell in the half-space occupied by the ions remains unchanged. In the 
second stage, at a sufficiently small density the distances between the molecules jump to 
a larger value (modification of the unit cell but not of the size of the molecules). Then the 
two stages repeat themselves until the value of the density is reached, where a transition 
from n-molecules to (n + l)-molecules occurs. After this transition the first two stages 
repeat themselves again, but now with ion configurations containing (n + l)-molecules, 
and so on. For a given potential U there is a largest admissible size N(U) of the molecules 
that can be formed. After the iV(£/)-molecules have been formed the remainder of the 
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sequence consists of the first two stages repeating themselves indefinitely. 

These phenomena were interpreted as consequences of the fact that for small U the 
n-molecules have only one bound electron, the other n — 1 electrons being diffusive. The 
g.s.c. appears then to result from the competition between a strictly convex effective 
repulsion between the molecules and the pressure exerted by the "free" electrons. 

On the other hand, for any density p, p c (U) < p < 1 — p c (U), the ground state 
configuration is periodic most homogeneous and the system is an insulator. Furthermore 
in the plane of the chemical potentials the density in the neutral domain appears to follow 
the Aubry sequence and the sizes of phase domains show the devil's staircase structure. 

We pass now to the non-neutral systems. In this case we have found that for p e < 



Pib(U), where b(U) is given by (|29|) the g.s.c. is the segregated configuration and the 
system is a conductor. The form of b(U) has been checked to be in agreement with our 
numerical results, the numerical results given in |J and the bounds given in [Tj| . 

For densities p e and pi such that Pib{U) < p e < 1/2, the set of periodic g.s.c. decom- 
poses mostly into sequences for which p% = np e , 1 < n < N(U). Given such a sequence, 
then for any pi between Pm in (U) and J oJ} iai ([/), the ground state consists of n-molecules 
which are most homogeneously distributed over the lattice, with one electron per molecule. 

For general values of (p e , pi) we came to the conclusion that for any electron density 
p e = p/q, with p relatively prime to q, there are periodic ground states with ion densities 
Pi = p/q, where p is any integer satisfying p' < p < p". If p« is such that p/q < Pi < 
(p+ l)/q, the ground state configuration is a mixture of the two phases with density p/q 
and (p + l)/q (for p = p" it is a mixture of the phase p"/q and the full configuration); for 
pi < p 1 jq it is a mixture of the phase p'/q and the empty configuration. In any case if the 
ground state is periodic, then the system is an insulator. 

The restricted phase diagram presents the following general properties. The non-trivial 
part of the (p e , pi) plane, i.e.— (U+2) < p e < 2, decomposes into three connected domains 
D_, D + and D, corresponding respectively to the empty, full and periodic translationally 
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non-invariant configurations. The domains D_ and D + are separated by domain D, 
contained in the vertical stripe [p* e ,—{U + /i*)], where p* e = p* e {U) < min{—2,—U/2}, 
and by two curves which correspond to the segregated phase (see Fig. 12). There are 
two partitions of D into connected subdomains that are of importance in the analysis of 
the restricted phase diagrams. In the first partition there appear subdomains D n , where 
the ions in the g.s.c. form n-molecules with one electron per molecule (pt = np e ). In 
each domain D n the ion density follows the Farey tree order (Aubry sequence) and the 
domains of distinct periodic configurations exhibit the devil's staircase structure. The 
second partition consists of the subdomains D Pe , where the electron density has a definite 
rational value p e = p/q. These domains form "curved stripes" going across D from 
the boundary with D_ to the boundary with D + . The ion densities in these domains 
constitute the sequence p'/q, (p' + l)/q, . . . ,p"/q, for some p' < p" that depend on p e 
and U. 

Finally, we have seen that theorem 1 of is partially incorrect and at the same time 
seems to be more general. We conjecture that this theorem should be valid for all U in 
the form: Given p e = p/q, with p relatively prime to q, and pi = p/q, then if the ground 
state is periodic, it is given by the Freericks-Falicov method. 
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